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Abstract. We investigate the evolution of the chiral magnetic instability in a protoneutron 
star and compute the resulting magnetic power and helicity spectra. The instability may 
act during the early cooling phase of the hot protoneutron star after supernova core collapse, 
where it can contribute to the buildup of magnetic helds of strength up to the order of 10^^ G. 
The maximal held strengths generated by this instability, however, depend considerably on 
the temperature of the protoneutron star, on density huctuations and turbulence spectrum of 
the medium. At the end of the hot cooling phase the magnetic held tends to be concentrated 
around the submillimeter to cm scale, where it is subject to slow resistive damping. 
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1 Introduction 

The origin of the magnetic held strengths observed in neutron stars and magnetars (highly 
magnetized neutron stars [1]) up to 10^^ G [2, 3] is still under debate. The most popular 
explanations include adiabatically compressed fossil helds of the parent star and dynamo 
generated helds [4]. More recently it has been suggested that the magnetic held of magnetars 
is related to a chiral asymmetry of particles, produced during the core collapse of supernovae 
[5]. An imbalance in the number of right- and left-handed fermions was previously studied 
in the context of QCD plasmas [6, 7] as well as applied to the early universe as a possibility 
to explain the generation and evolution of cosmological magnetic helds [8, 9]. This so-called 
chiral magnetic effect or chiral magnetic instability was also suggested to account for the 
observed kicks that accelerate neutron stars [10]. 

In a pure electron-positron plasma, the chiral magnetic instability does not allow the 
growth of seed magnetic helds [11], but in the presence of neutrinos an electroweak plasma 
with neutrino-antineutrino asymmetries was found to be able to amplify magnetic helds to 
interesting scales for neutron stars [12, 13]. Another crucial ingredient to take into account 
are the spin hip interactions due to the hnite electron mass which violates chirality. This 
tends to decrease the asymmetry between left- and right-handed electrons faster than it is 
created by electroweak processes [14]. In addition, it was claimed that the chiral asymmetry 
in the forward scattering amplitude of electrons off nuclei due to the electroweak interaction 
can create a magnetic held instability in the same way that the chiral asymmetry does, but 
which acts on much longer time scales and is not washed out by chirality-hipping processes 
[15, 16]. 

It is not yet clear if the chiral magnetic instability can transfer sufficient energy stored 
in chiral fermions into magnetic held energy to give a signihcant contribution to the magnetic 
helds inferred for neutron stars and magnetars. In the present paper we aim to model the 
chiral magnetic effect in this environment, review the underlying assumptions proposed for 
this mechanism to work and understand its physical implications. We solve the evolution 
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equations for the chiral chemical potential, the chemical potential of the background species 
on which electrons scatter, the magnetic energy and the magnetic helicity power spectra. This 
allows us to estimate the conditions for which magnetic fields can be amplified in a neutron 
star. We show that a seed magnetic field can be amplified to small scales shortly after the 
collapse in the core of hot stars through the chiral magnetic effect and that for its surface or 
cooled down neutron stars this mechanism is not effective to generate the strong magnetic 
fields observed. 

The remainder of the paper is structured as follows: In Sect. 2 we summarize the condi¬ 
tions in a protoneutron star relevant for magnetic field evolution and set up the basic modified 
MHD equations. In Sect. 3 we solve the evolution equations, estimate the maximal magnetic 
field strength and discuss assumptions and uncertainties. We summarize our results and con¬ 
clude in Sect. 4. Throughout the paper we will use Gaussian natural units, cq = ^ = /cs = 1, 
and the electric permittivity and magnetic permeability of the vacuum are set to cq = 1/ (dvr) 
and Ho = 47r, respectively. 

2 Framework and basic equations 

2.1 Thermodynamics of a protoneutron star 

Immediately after core collapse, a protoneutron star reaches temperatures of the order of tens 
of MeV in its core. To a given temperature corresponds a chemical potential = Hn— fJ'p = 
He — Hv [17]) that can be used together with the fact that neutrinos are trapped inside the 
core at this stage, such that the lepton fraction Yl is temporarily conserved. The relation 
Flu-b = Ue + Uu, where ub = Tin + np \s the baryon number density, and electric neutrality. 
Tie = Up, allows us to estimate the number densities and chemical potentials of the particle 
species involved. 

When a massive star collapses protons are converted into neutrons by capturing left- 
handed electrons cl + p ^ n + Ve^ i producing an asymmetry between the number of left- 
and right-handed electrons = Nl — Nb- Such electroweak reactions are known as URCA 
processes and their emissivity is [18] 

457yr 

euRCA = ecGpmnmpHeT^ , (2.1) 

where gA — 1-26 is the axial-vector coupling of the nucleon, 9c ~ 0.24 is the Cabbibo angle, 
Gp = 1.166 X 10“® GeV“^, and irip are the masses of the neutron and proton, respectively 
and He is the electron chemical potential. Since He T, the rate of electron capture is then 

P _ euRCA 
^ ~ HeYLUB 

If the URGA processes are not in thermodynamic equilibrium with the inverse reactions, an 
asymmetry can build up. This is the case, for example, if neutrinos escape the neutron 
star, which occurs when their mean free path is larger than the neutron star radius. This 
condition is met beyond the neutrino sphere or when enough time has passed for the star to 
cool down to the point when it becomes transparent to neutrinos, roughly 10 seconds after 
collapse. 

For the typical momenta of the particles in the core of a protoneutron star electron 
capture takes place, while in the crust or when proton and electron concentrations are low 
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momentum conservation highly suppresses electron capture and an additional particle is re¬ 
quired to absorb momentum, such as another proton or neutron. If pF,n > PF,e + PF,p-, the 
previous rate is modified to [19] 
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where ~ 15 is the pion-nucleon fine structure constant. 

During the hot initial phase, the charge carriers in the neutron star are semi-degenerate. 
There are no simple equations for the conductivity for this case. In the degenerate limit the 
conductivity is given by the following expression [20], 
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with Pp the proton density. In the relativistic non-degenerate high temperature limit the con¬ 
ductivity would be dominated by the pair plasma and be of the order of the temperature [21] 
which is lower than the obtained through (2.4). However, it turns out that the final magnetic 
field is insensitive to the conductivity in this range of values because the magnetic field grows 
on timescales much shorter than the dynamical timescale of the system and then saturates. 
Therefore, for our calculations we will use eq. (2.4) for the conductivity, which is closer to the 
relevant conditions in the protoneutron star, and we will also assume the temperature T and 
the conductivity a to be constant. 

Even though in the core of a collapsing supernova electrons are relativistic, the fact 
that they are massive suggests that we should not take them as strictly chiral particles, 
since the amplitude of a positive helicity component for a left-chiral state is approximately 
{E + me — p)/{E + me) — {me/E). This means that there is a probability ime/E)^ that 
a scattering electron of a certain chirality flips into the opposite chirality state - either by 
Rutherford scattering, electron-electron scattering or Compton scattering - which tends to 
decrease Rutherford scattering dominates in this case, which allows us to write the 
chirality-flipping rate as {E T) [14] 

( 2 . 

where e is the electron charge and me is the electron mass. 
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2.2 Evolution equations 

When a chiral imbalance is present, such as the one originated by electron capture, the Adler- 
Bell-Jackiw anomaly implies that a current 

J5 = = Jl - Jr (2-6) 

will be induced, whose partial derivative, instead of vanishing, is related to the Chern-Simons 
current through 

« = = dpK >^. (2.7) 

With = J d^r'075?/; and Ncs = f d^rAT*^ space integration implies the conservation relation 

— (iVs - Ncs) = 0 . 
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Thus the Chern-Simons number of the electromagnetic held connects the chiral asymmetry 
to the magnetic helicity 


dt 
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with the magnetic held B and vector potential A. 

Maxwell’s equations are then modihed by the introduction of a current contribution 
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in the presence of a chiral imbalance, with chiral chemical potential /is = (/i£, — fiji)/2. This 
affects the magnetohydrodynamics (MHD) equation that now takes the form 


dtB = V X (v X B) + r/AB 


vr 


T/ZisV X B, 
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with rj = l/(47rcr) being the resistivity. In the following we will neglect the velocity held v 
so that our subsequent analysis applies in the plasma rest frame, provided that the velocity 
held is sufficiently smooth, which we elaborate in appendix A. It will be left to future work to 
investigate under which conditions this is a good approximation in the presence of turbulence 
and other contributions to the velocity such as rotation. 

In Fourier space and introducing the expansion of the magnetic held into a left- and 
right-handed part 


B(k) = 6+h+ + i,^h. 
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where e is an arbitrary unit vector perpendicular to k, one can then rewrite (2.11) as 


dtb^ = -rjk ( k ± • 


(2.13) 


The magnetic held energy density and helicity density can be written in terms of the 
power spectra in Fourier space, = A:^|B(k)p/2 and = —Airiklk. x B(k)].B*(k) , as 
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Combining (2.12) and (2.13) and multiplying with the magnetic held complex conjugate, the 
power spectra evolution is given by 


dtPm = -^ J din kk‘^ ^2Mk + , (2.15) 

dth = —^ J din k{2k‘^Hk + 32e^p,5Mk), (2.16) 

where V is the volume. We can translate particle number into chemical potential using 

^5 = {pI + , (2-17) 
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which can be approximated to hnear order in /rs to = c{T, with 

m2 7^2 

c(r,Me) = ^ + ^. (2.18) 

TT^ 6 

To express the evoiution of the chirai chemicai potentiai we take into account the processes 
described by (2.5) and (2.9) that affect the number of ieft- and right-handed particies, resuiting 
in 

e 2 

- 2Tf , (2.19) 

where is the equilibrium value of /rs in the absence of resistivity. This term represents an 
effective chemical potential generated by the interactions of electrons with background species 
such as neutrinos that act as sources of the asymmetry, thus containing the term proportional 
to (2.2). We can roughly estimate it by considering the number Nb of background particles, 
such that the processes that change N^, neglecting the magnetic held contribution for now, 
can be written as 


= ±r^Nb - 2r;iV5 • 

Comparing this with (2.19), results in 

( 2 . 20 ) 

Nb = 21/c(r,/ie)^|/i5,fe| • 

-L w 

( 2 . 21 ) 


Furthermore, we approximate the background particles at the temperatures of the core of a 
protoneutron star to be non-degenerate relativistic fermions with Qb degrees of freedom. This 
is plausible since the chiral asymmetry mostly results from URCA processes involving chiral 
neutrinos whose chemical potential is at most of the order of the temperature which would 
give rise to order one corrections. From thermodynamics we can then relate the number of 
background particles with the temperature using Nb = 3C{3)VgbT^/, yielding 


_ ^ F^ T3 
87r2 ^^F/c(T,/ie) ■ 


( 2 . 22 ) 


This effective background chemical potential can also include a possible contribution from the 
difference of the forward scattering amplitudes of left- and right-handed electrons on nucleons, 
that was considered in the form of an effective potential V 5 in [15] and [16]. Requiring energy 
conservation of the combined system consisting of the magnetic held, the chiral asymmetry 
in the electro sector and the background particles implies an evolution equation for g^^b of 
the form of 

dtg'5,b = , (2.23) 

^^5,b 


as we will see below. 


2.3 Magnetic field amplification 

The initial magnetic held of a protoneutron star, which for example can result from adiabatic 
compression of the stellar seed held during collapse, is affected by the chiral magnetic effect 
in a way which depends on the scales we are interested in. From (2.13) for growing modes we 
obtain the condition 

2 e 2 
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such that magnetic field modes with k > decay due to resisitivity. For k < k^ the magnetic 
field mode with the same sign for helicity and //s also decays while the mode with helicity 
signs opposite to grow. All magnetic field modes are damped due to finite resistivity 
with the resistive damping rate F,. = rjk'^ whereas the growth/decay rate due to the chiral 
instability is given by 

r^(A:) = —r//c|/r5| = . (2.25) 

vr K 

The maximal total growth rate Ftot = F^ — F^ occurs at Fmax = which corresponds 

to the wavenumber ^ 5 / 2 . 

Let us analytically analyze the expected behavior of /is by setting = 0 in (2.19). If we 
normalize to the maximal value that the helicity can take, FFmax(fc) = SnMf^/k, this gives 


Fj/is^fc - 2r/e^[7rl/c(T,/ie)] V dlnfcA:Mfc(FFfc/FFmax) 
Tf + 4r/e4[7r2c(T, /ig)]" Vm 


(2.26) 


When the magnetic field is negligible, ~ /is^s, and modes smaller than k^ grow exponentially 
at the rate Ftot- The magnetic field terms begin to dominate when pm ^ 7 r^c(T, /ig)/(4r/e^)r j. 
In this limit the flipping rate is negligible compared to the magnetic field induced rate and 
the instability yields 
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In general 7 ^ /rs s and the terms of (2.19) compensate each other. As a consequence the 
magnetic helicity density will change linearly with time at a rate 


^ . (2.28) 

The fact that we have almost maximal helicity implies that also the magnetic energy density 
changes linearly with time, either growing or decreasing according to the sign of (p^ — P 5 ^b)/h. 
The helicity would be constant only if Fj = 0 or if /is = /isAssuming that the magnetic 
energy is concentrated around a characteristic scale ko = kb{p^) and that the helicity is max¬ 
imal and has opposite sign to /is implies that also pm will be constant. In this case the chiral 
magnetic instability reaches saturation, where the growth and damping rates compensate 
each other. 


2.4 Energy balance 

By definition of a chemical potential the energy associated with the chiral asymmetry is 
given by dEs = pbdN^ and with Es = 0 for /is = 0 results in 



c{T,Pe) 2 


(2.29) 


Differentiating this with respect to time, using (2.19) with the helicity normalized to its 
maximal value and inserting k^, we obtain 


dtP5 = - 2 t / [ dlnk^ (hksgn{p5) 

. V \ iimax 


+ H] ~ 2c(r, /ie)r//is(/r5 - 1^5,b) • (2.30) 


We can also estimate the change in magnetic energy that a finite /is can induce by using that 
the instability produces maximally helical fields and that, as we have seen, the growth has 
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its peak at /C 5 / 2 . Thus dEm — A: 5 |d^|/( 87 r) and using (2.9) gives dEm — Vc{T, fj,e)fi 5 dfj, 5 . If 
there is an initial chiral asymmetry the increase in magnetic energy density is given by 

52^(4-m1). (2.31) 

The total energy density 

Act = P5 + Apm , (2.32) 

then only depends on the initial value /is* which implies that the maximal increase in magnetic 
energy density obeys Apm < Ptot- 

Eq. (2.15) can be rearranged into 


dtPm = -^rj j d\nkk^^ + ^sgn(^5) 


Hk 


which together with (2.30) provides the rate of change of the total energy 
dtPtot = dtPm + dtP5 


= —2p f dlnk^^ 


(k - k^f + 2k^k ( 1 + sgn(^5) 
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(2.34) 


The term proportional to is responsible for the energy exchange with external particles 
and we see that apart from it, since the integrand of (2.34) is non-negative, the energy 
decreases due to chirality flips and hnite resistivity. It is strictly conserved only for the case 
P5 = P 5 ,bi and if the magnetic energy is concentrated in the mode k^ and helicity is maximal 
with opposite sign to which are the same conditions mentioned above for dtPm = 0 . 

Similarly to the energy associated with the chiral particles, the energy density pb asso¬ 
ciated with the background species is given by dEb = p^^bdNb- Using (2.21), this gives 


Pb = 


Eb 

V 



(2.35) 


The initial chemical potential of the background species, pb^bi, gives us a measure of the 
maximal energy that can be transferred into magnetic energy density 


^Pm < c(r, Pe)^plbi ■ (2-36) 

For the interactions to conserve energy one has to set dtPb = —dtPh in the absence of magnetic 
helds and, computing the time derivative of (2.35) and using (2.30), this yields the evolution 
equation for the background species (2.23). This shows that p^^b changes typically with the 
rate T^ and that will be in equilibrium for p^ = p^^b when the magnetic held is concentrated 
around k^{p^^b) nnd has maximal helicity of sign opposite to p^p, given by (2.22). Fig. 1 
shows the dependence of the ratio T^/T f on the temperature for a protoneutron star. 

Once the energy contribution (2.36) is added to (2.34), the only remaining source of 
energy change is resistive damping. Inserting (2.22) into (2.36), we obtain 
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max 
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3C(3) 

87r2 ^fc{T,pe) 


~ 8.3 X 10“^ 
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(2.37) 
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Figure 1. Left panel: The ratio T^j/r/ obtained using (2.2) and (2.5). Right panel: Estimate 
maximum magnetic field amplification i?niax = \/due to the chiral magnetic instability as a 
function of temperature. Logarithmic units to base 10. 


where in the last expression the degrees of freedom of the background species were taken 
as <76 = 2. (2.37) can be used to predict the maximal held strength generated for a given 
temperature and is shown in hg. 1. We can see that the maximal magnetic held amplihcation 
strongly increases with temperature. 

Up to now we assumed that 7 x 5 evolves only in time, neglecting the possible contribution 
of its spatial evolution, as studied by Ref. [22]. There it was shown that the system can 
be unstable with respect to growing inhomogeneous modes of 7 x 5 if fe^/4 > 3 e^i?^/( 87 r^T^) 
[22]. If we consider the fastest growing mode ^ 5/2 and the maximum magnetic held strength 
generated from (2.36), we can estimate whether the inhomogeneity of /is plays a role or not 
in our case. We hnd then that the solution of (2.19) will be stable if r^o/T f < d/vr, which is 
verihed for the temperature values we are interested in, as is clear from the left-hand side of 
Fig. 1. This implies that for our purposes = /i 5 (t) is a safe and justihed assumption. 


2.5 Density fluctuations 

Neutrinos of energy Ei, ~ 3T are trapped at the high temperatures in the protoneutron star 
core, where densities easily reach ub = 2no, with no — 1-7 x 10^® cm“^ the nuclear matter 
number density. The mean free path for absorption by a neutron is [23] 
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and the mean free path for scattering with a neutron is 
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(2.39) 


The absorption mean free path is more important than scattering and yields ~ 1.5 m for 
T = 20 MeV and ~ 10 cm for T = 40 MeV. The scattering mean free path yields ~ 1 m at 
20 MeV and ~ 14 cm at 40 MeV. These typical temperatures will be used in the following 
section. 
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To take into account the fact that the interior of a young neutron star is turbulent we 
consider the existence of density fluctuations 5p relative to the average density p. These 
density perturbations seem to amount to at least 25% [24], If the scale of the fluctuations is 
smaller than the neutrino mean free path, locally these regions can amplify a seed magnetic 
field, since the neutrinos stream freely on that scale such that the URCA processes and 
their reverse processes are not in thermal equilibrium and the rate of production of chiral 
imbalance will be of the order of the direct URCA rates T^. To study the influence of the 
fluctuations somewhat more quantitatively we introduce an effective creation rate of chiral 
imbalance T^ = F^dp/p which is a rough estimate of the difference of the absorption and 
emission rates of left chiral electrons due to the electroweak URCA interactions. Rewriting the 
chiral asymmetry equilibrium value (2.22) and the characteristic wavenumber of the instability 
(2.24) in terms of the density fluctuations, we have 


6pr^ ,_47MT^f6pr^.Y^ 

^5,6 87r2 T/c(r,;Ue)’ 2>e^Q{‘i)gb \pFf ) 


(2.40) 


From (2.37), the resulting maximal field amplification with respect to the density fluctuations 
then becomes 
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(2.41) 


3 Solutions of the evolution equations 


We now apply the previous treatment to the core and the neutrino sphere of a protoneutron 
star. In a core collapse supernova with a progenitor mass ~ 8Mq, as described by [25], shortly 
after core collapse the lepton fraction is Yl ~ 0.3. For the chemical potential difference we 
consider the two realistic cases Ap = 80 MeV and 60 MeV, which correspond to temperatures 
of ~ 40 MeV and 20 MeV, respectively, for a core density of 2no, as before. We can then 
compute the number density and chemical potential of each species as described in Sect. 2.1 
which results in the electron chemical potential Pf, ~ 260 MeV and the proton densities 
Pp ~ 1 . 3 x 10 ^^ gcm“3 for 40 MeV and 1 . 2 x 10 ^^ gcm“3 for 20 MeV. From this the conductivity 
is obtained from eq. (2.4). 

We solve the system of eq. (2.15), (2.16), (2.19) and (2.23) for 90 wavenumber modes 
spanning from femax = 2 A :5 to /Cmin = lO^'^fcs, with constant width in logro time scale 

used is normalized to the resistive damping time of the instability 


idamp — 


32TT^c{T,pe) / dpT^A ^ 
%{3fe‘^glp \p Ff ) 


(3.1) 


We first analyze the case for which the density fluctuations are close to the average density 
~ p and then consider that 6p\s lower than p by one order of magnitude. After the initial 
10 seconds of its life, the neutron star becomes transparent to neutrinos and the magnetic 
field amplification can be estimated by putting 8p/p = 1. 

In fig. 2 the evolution of the chiral magnetic instability with a vanishing initial value for 
the chiral chemical potential is shown and table 1 contains the respective |/r 5 ,;)|, F”^, FJ^, 
idamp and values taking 5p/p = 1. 

From (2.41), the maximum field strength that can be reached by the instability is inde¬ 
pendent of the initial magnetic seed field. Its strong temperature dependence reflects the fact 
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Figure 2. Time evolution of the chiral chemical potential normalized to the equilibrium value, 
M5/IM5 ,&|j relative difference of the chiral chemical potential to the equilibrium value, (/is —/i5,6)/|/i5,b| 
and, in logarithmic units, relative deviation of the helicity density from its maximal and minimal 
value, 1 ± /i//imax- The left panel is for a temperature of T = 40 MeV and seed field Bq = 10^^ G, 
and the right panel is for T = 20 MeV and a seed field of Bq = 10® G. 



Figure 3. Time evolution of the magnetic energy density Pm and total energy density ptot- Also 
shown is the initial total energy density which limits the maximal magnetic energy density that can 
be reached by the instability. In the left panel T = 40 MeV and in the right panel T = 20 MeV. 


T (MeV) 

M (MeV) 

T)))' (s) 

r/' (s) 

^damp (s) 

h ^ (cm) 

Smax (G) 

40 

2 X 10-3 

1 X 10"® 

3 X 10-32 

6 X 10-3 

3 X 10-® 

8 X 1034 

20 

4 X 10“® 

9 X 10-3 

3 X 10-32 

0.04 

1 X 10-3 

1 X 1033 


Table 1. Equilibrium chiral asymmetry (2.40), chiral asymmetry creation and depletion rates 
(2.2) and (2.5), respectively, damping time (3.1) and characteristic scale (2.40) values for Sp/p = 1. 
Maximal magnetic field amplification computed using (2.41). 


that Flo has a stronger dependence on temperature than Tf , which means that the higher 
the temperature, the higher the magnitude of ^5 and the sooner the chiral magnetic insta¬ 
bility develops and subsequently damps. The values predicted for the maximal field strength 
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(2.41) are listed in the last column of table 1. For the examples we chose magnetic seed helds 
Bq = 10^^ G for T = 40 MeV and Bq = 10® G for T = 20 MeV to illustrate how the magnetic 
held amplihcation is limited by energy conservation when the magnetic energy gets close to 
the maximum allowed value. In hg. 3, for T = 40 MeV the magnetic energy density grows 
steeply up to close to the total energy of the system, corresponding to i?max — lx 10^^ G, 
within a few ^s and subsequently the magnetic held decreases by a factor of more than 10 
within a few seconds. For T = 20 MeV, the magnetic held grows within about 4 seconds up 
to Smax ^ 1 X 10^2 G. 

The total energy density 


Ptot — g {pb Ptii 


(3.2) 


includes the energy density from (2.29), the energy corresponding to the background 
particles coupling to the chiral electrons from (2.35) which ensures that the total energy 
due to the scattering terms is conserved, and the magnetic energy density pm- The initial 
total value is not exceeded and, as predicted, ptot decreases only due to resistive damping 
following the dissipation of pm- The number of wavenumber modes considered is sufficiently 
large for the decay of the magnetic energy to correspond to a smooth curve. The decay 
changes pm roughly linearly with time, as expected from the discussion in Sect. 2.3. 

The chiral asymmetry is built up through the capture of left-handed electrons until an 
equilibrium with the spin-flip processes is reached at ps = When the magnetic held 

starts to be amplihed, the term in (2.19) proportional to the magnetic helicity will eventually 
dominate. At this point the asymmetry ps will start to decrease as chiral energy is transferred 
into magnetic energy. 

The curve (ps — Pb,b)/\Pb,b\ in hg- 2 is very close to zero after the equilibrium value is 
reached, which implies that when the magnetic held terms begin to dominate the evolution of 
Pb and occurs in lockstep. The chiral chemical potential grows and reaches equilibrium 
at a value close to until the magnetic held term starts to dominate and depletes the chiral 
asymmetry. The scattering of electrons, covered by the term proportional to p^^h does not 
allow for p5 to be replenished. 

The magnetic helicity density, normalized to the maximal value /imax(fc) = {87r/V){Mk/k), 
depends on the mode considered and maintains its initial value, here simply chosen as 
h® = hmax/2, until the amplihcation of the magnetic held makes it either grow to its maxi¬ 
mum or decay, if the sign between helicity and p^ is the opposite or equal, respectively. Thus, 
helds amplihed by the chiral instability turn into maximally helical helds. Furthermore, once 
magnetic held growth sets in the evolution is essentially independent of the initial helicity. 

It is also interesting to compute how the magnetic held power spectrum evolves with 
time. Fig. 4 shows the time evolution of the magnetic held power spectrum for a hat and 
a Kolmogorov initial spectrum. As expected, the hnal magnetic held power spectrum is 
not very sensitive to the initial magnetic held power spectrum. The magnetic held power 
spectrum peaks at wavenumbers close to k^/2, while it decays with time due to resistive 
damping for k > kb- Since k^ is proportional to the evolving chiral chemical potential p^, 
see (2.24), which decreases for log;^o(V^damp) > 2, hrst steeply and then smoothly, see hg. 2, 
with growing time the peak in the magnetic power spectrum moves to smaller k. This is 
rehected in hg. 4 which also shows that the total magnetic energy grows exponentially for 
times ~ lO^tdamp ^ ^ ^ lO^tdamp) then saturates and gets damped for subsequent times. 
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T = 40 MeV, Bo = q 



Figure 4. Time evolution of the magnetic field power spectrum normalized to the initial magnetic 
energy density, as a function of wavenumber k normalized to k^. The power spectra are 

shown for equally spaced intervals in the logarithm of time between t = tdamp and t = 10®tdamp) for 
T = 40 MeV. Left panel: Initially fiat power spectrum. Right panel: Initial power spectrum has a 
Kolmogorov distribution. 


In the simulations of the magnetic field power spectra we considered the most relevant 
modes k for magnetic amplification: The peak of the magnetic field power spectrum yields 
a maximal growth for k^/2 and taking into account much smaller wavenumber modes up to 
the size of the neutron star radius, k ~ 1/(10 km), does not significantly change our results. 

We can also estimate the time dependence of k^ in the damping regime: Eq. (2.15) shows 
that amplification stops and resistive damping sets in when 2r]k‘^t ~ 1. Therefore, we expect 
the scaling 

^5 ~ ^5 ( ^ 1 ) A^5 ~ /^5 ( ^ I ) (3-3) 

where k^ and are the values of k^ and /is, respectively, at t ~ tdamp- Fig. 4 allows us to 

estimate the length scale at which the power spectrum peaks after the magnetic field growth 
ends. For T = 40 MeV we obtain k~^ ~ 0.1 mm while for T = 20 MeV we find k~^ ~ 3 

cm. This is significantly smaller than the neutrino mean free path discussed in Sect. 2.5. 

Therefore, in the presence of significant density fluctuations on these length scales we expect 
that magnetic field growth due to the chiral magnetic instability is possible. 

Let us finally turn our attention to the case involving density fluctuations dp/p = 0.1 
in the core of the neutron star, as illustrated in fig. 5 for seed magnetic fields of 10® G. For 
T = 40 MeV one obtains a field amplification of now only ~ 1 x 10^^ G, whereas for 20 MeV 
it yields ~ 2 x 10^^ G. 

3.1 Neutrino sphere and cold neutron star 

If we consider the typical radius of a neutron star to be 10 km, the density at the neutrino 
sphere is ~ 10^^gcm“^. The lepton fraction can be roughly taken as ~ 0.1 [26] and 
the average neutrino energy is at most 16 MeV [27], which translates to a temperature of 
T ~ F^i//3 ~ 5 MeV. At this and lower temperatures, an additional particle is required for 
electron capture to occur, similarly to the modified URGA process NYp + e^ —?• N + n + zze^, 
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Figure 5. Energy densities obtained when the factor 5p/p = 0.1 is included, accounting for the 
density fluctuations in the neutron star core. 


where N can be either a neutron or a proton. The estimate (2.22) gives a chiral asymmetry 
of 1/^5,fel 10 MeV which in turn through (2.41) gives negligible maximal magnetic fields. 

Let us now briefly concentrate on mature neutron stars of about 10® years old, where 
the core temperature drops to typically 2 x 10® K and the main cooling mechanism is surface 
photon emission [25]. The URCA rate yields an also very small instability equilibrium value 
\^J‘ 5 ,b\ ~ 10“®^ MeV, which in turn corresponds to a field growth only around extremely small 
k and in extremely large time scales, indicating that the chiral asymmetry is also not effective 
in generating magnetic fields for this regime. We, therefore, obtain significant magnetic field 
amplification only during the hot initial phase of the neutron star. This is in contrast to 
[15, 16] which included the difference in the forward scattering amplitudes for left and right 
chiral electrons which gives rise to a constant potential term V5. These authors effectively 
introduced the rescaling —)• /rs + V5 in the terms involving the magnetic field but not 
in the terms for the chiral asymmetry evolution due to interactions with the background, 
yielding an artificial steady source of magnetic energy. However, in our opinion this rescaling 
should be performed in all terms so that it only introduces a constant shift of which can 
be eliminated by redefining p^. In other words, the thermodynamic equilibrium abundances 
can only depend on the total chiral energy difference P 5 + V 5 , and not only on p^. As a 
result, Vsshould affect neither the scale nor magnitude of the magnetic field amplification. In 
contrast, if only the kinetic part of p 5 would be given by (2.22), substituting p^^u —>■ p^^u + Vb 
in (2.36) and using p^^i^ ~ /is ^ with (2.22) would yield 


Pm < c(r, pe) 




rp3 

ciT, Pe) 



(3.4) 


But since T^ oc G\ and oc Gp this would not vanish in the limit of Gj? —)• 0, as is expected 
since if parity is conserved and there is no chiral asymmetry, the energy associated with the 
chiral asymmetry should vanish. This indicates that the evolution of the system should not 
depend on V5. The consideration of a term of this kind has been also discussed and discarded 
in [14, 29] under similar reasoning. In Ref. [28] the authors claim that the maximal magnetic 
field energy density Req/(^ tt) is given by the thermal energy of the nucleons and electrons but 
no derivation is given and it is unclear how this is related to V 5 . Furthermore, saturation of 
magnetic field growth is introduced ad hoc by substituting + V5 —?• {p^ + V5 )/(l + 
without derivation. 
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4 Summary and conclusions 


In a supernova core collapse electron capture creates an imbalance between left- and right- 
handed electrons /rs ~ eV - keV. It has been suggested that the chiral anomaly can transform 
the energy associated with this chiral imbalance into the growth of helical magnetic helds, 
possibly up to the high values that have been observationally inferred for neutron stars and 
magnetars. In the present work we investigated this possibility within a semi-analytical 
approach with specihc emphasis on the evolution of the total energy which can only decrease or 
stay constant and thus limits the maximally possible magnetic held strength. While neutrinos 
are trapped in the core, density huctuations allow for local thermodynamic disequilibrium 
between URCA and inverse URCA rates due to neutrino free-streaming on sufficiently small 
length scales ~ cm, which prevents this imbalance to be washed out by the inverse reactions. 
For length scales 10“®cm < 7r/(2e^|/r5|) ^ i,/ ^ 15cm the chiral magnetic effect can then 
create magnetic helds of roughly maximal helicity on time scales short compared to the 
evolution of the neutron star before they saturate due to the limited energy associated with 
the chiral lepton asymmetry. 

For a core temperature of 40 MeV, we obtain maximal magnetic held strengths Umax ~ 

G on tens of nanometer length scales reached within microseconds. For lower temper¬ 
atures, such as 20 MeV, the magnetic helds are smaller and concentrated on larger length 
scales, and the growth rates are lower. This suggests that the range of held strengths and 
power spectra due to the chiral magnetic instability depend on the initial temperature of 
the protoneutron star. The generated helds are not strong enough to account for typical 
magnetar held strengths and tend to be produced on submillimeter length scales rather than 
dipolar helds on the linear scale of the star. We also hnd that outside the neutrino sphere, as 
well as in a cold neutron star at temperatures below ~ 10 MeV, the chiral instability can not 
lead to signihcant held amplihcation, due to the fact that lower temperatures imply smaller 
asymmetry values 

We briehy summarize the differences between our study and other recent work on the 
chiral magnetic instability in neutron stars. Maximum surface helds of 10^® G were estimated 
in Ref. [5] by considering a very high and constant /U5 = 200 MeV. The procedure in Ref. [14] 
yields a chiral asymmetry ~ 10“^^ MeV (for T = 30 MeV), several orders of magnitude 
lower than our result motivated by the electroweak electron capture rate. The responsible 
mechanism for the magnetic held growth in Ref. [15, 16, 28] is stated as being due to a potential 
term V5 that accounts for the parity asymmetric forward scattering of chiral electrons and 
nucleons and which acts on a much longer time scale, being relevant for cold neutron stars. 
In our treatment any chiral asymmetry in the forward scattering of electrons on background 
species does not separately contribute to the enhancement of the magnetic held because only 
the total asymmetry energy /xs -|- V5 should enter the evolution equations, as also concluded 
in Ref. [14, 29]. We rather believe that the magnetic held evolution only depends on the 
asymmetry between left- and right-handed electron abundances which in turn is a function 
of the ratio of electroweak URGA and spin hip rates. 

One important approximation taken throughout this work was neglecting the role that 
turbulence may play in suppressing the chiral instability by setting v = 0 in the MHD 
equation. This should be a good approximation as long as the velocity held is sufficiently 
smooth on the instability length scales or the power index sufficiently large, as shown in 
A, which allows to transform into an inertial frame moving along with the plasma. We 
also consider the temperature and the resistivity of a protoneutron star to be constant over 
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the first initial stage of evolution after the supernova collapse. The later should be a good 
approximation since for T > 20 MeV all relevant time scales, including the instability growth 
scale, are short compared to the scale on which the temperature changes, which is a few 
seconds [17]. Lower temperatures imply instability time scales longer than a few seconds on 
which the cooling of the protoneutron star should be taken into account. The resistivity that 
we employed in this work assumes that particles in the protoneutron star core are degenerate 
whereas a semi-degenerate regime is more realistic. This affects the resistive damping rate 
and thus the damping time and the timescale in which the instability grows. Additionally, 
since the instability timescale is in any case small compared to the dynamical timescale, the 
final state, in particular the finite magnetic field strength, is not influenced by this uncertainty 
because it is determined by saturation of the magnetic field energy at a value comparable to 
the energy in the chiral asymmetry and given by eq. (2.37). 
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A The role of turbulence 


In this section we show under which conditions the assumption that the fluid turbulence can 
be neglected is a good approximation. Comparing the first and third terms of (2.11) provides 
us with an estimate of how large the velocity in the core of a protoneutron star has to be 
to dominate over the chiral term in the MHD equation studied in the present work. Let us 
assume a fluid velocity spectrum of the form 


{v\T,k)) 




k 


(AT) 


with ki being the inertial wavenumber. In terms of the length scale i = 21^jk and of the 
root-mean-square velocity Uj-ms = \/ ('*^^(T, A:)), this gives us for the velocity flow 


= ■*^rms 



(A.2) 


with L the integral length scale and n the power index. From the MHD equation with the 
chiral anomaly (2.11), the velocity and anomalous term can be estimated as 


_ , vB 

V X (n X B) ~ —, 

V7 T3 

27 ^ 2 ^ ^5^ 27rVr 

such that the relative importance of the first is roughly dictated by 


(A.3) 


V X (n X B) 
e 2 /( 27 r 2 a)/r 5 V x B 


2aLvr 




-(n/2+1) 


(A.4) 
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by considering that the relevant scale for the instability is £ = 27r//c5 = (7r/e)^|/U5|“^. 

We take the example of a protoneutron star with T = 40 MeV, with a corresponding 
conductivity of a ~ 0.21 GeV computed from (2.4), that according to table 1 has an electron 
chiral chemical potential of 2 x 10“^ MeV. Considering the scale at which most energy will be 
concentrated after the hot cooling phase being L ~ km, for a Kolmogorov velocity spectrum 
(n = 2/3), we obtain 


V X (rr X B) 
e2/(27rV)^5V X B 


~ 5 X lO^Urms- 


(A.5) 


This implies that for n = 2/3, the effect of the chiral instability dominates unless the fluid 
velocity of the protoneutron star core is Urms ^ 10“^. As an upper limit to the typical 
velocities implied we consider the example of 4 x 10® cm/s [24], rendering Urms ~ lO”^^ which 
indicates that in this regime our results do not apply for Kolmogorov turbulence, but only 
for larger power indices, such as n = 4/3 in this case. For a lower limit on the fluid velocity 
we take the example of 10® cm/s [30], giving Urms — 3 x 10“®, which according to (A.5) 
shows that the MHD equation will be dominated by the chiral anomaly for a Kolmogorov 
type spectrum in this case. 

The velocity term in the MHD equation becomes less important as n increases: for a 
Kraichnan type spectrum (n = 1), (A.4) is of the order unity and smaller for Urms ^ 10“®. 
The previous estimates make clear that turbulence can be neglected for relatively small fluid 
velocities or large power indices, but that depends on the details of the velocity spectrum at 
play in the core of the protoneutron star. 
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